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The  AASERT  Grant  supported  two  students  at  Rutgers  University,  Mr.  Dov  Chelst 
and  Mr.  Kevin  Rosema. 

Dov  Chelst  received  his  Ph.D.  in  May  1999.  The  title  of  his  thesis  was  “Modified 
Two  Component  Plasmas  and  Generalizations  of  Schwarz’s  Lemma”  .  Dov’s  thesis  dealt 
with  statistical  mechanics  of  one  dimensional  systems  with  Coulomb  interactions.  This  is 
a  problem  on  which  some  beautiful  analytic  work  was  done  a  long  time  ago  by  Baxter, 
Edwards,  Lenard,  and  others.  Their  solutions  employed  an  implicit  equation  involving  a 
continued  fraction.  In  addition,  Lenard  used  the  analytic  theory  of  continued  fractions  to 
prove  the  existence  of  an  impUcit  solution.  However,  he  did  not  realize  that  this  theory 
could  also  be  used  to  prove  the  analyticity  of  this  solution  with  the  aid  of  an  implicit 
function  theorem. 

Chelst  thesis  extends  Lenard’s  analysis  to  include  a  large  class  of  systems  so  called 
modified  two  component  plasmas.  This  includes  both  Coulombic  and  non-Coulombic  types 
of  interaction. 

The  motivation  for  this  arose  from  some  experimental,  numerical  and  approximate 
anal3dical  calculations  in  higher  dimension.  These  indicate  that  Coulomb  systems  with 
hard  cores  undergo  a  rather  peculiar  kind  of  liquid-vapor  transition.  It  seemed  that  it 
would  be  useful  if  one  could  investigate  the  simpler  one  dimensional  system  exactly.  This 
Dov  has  done  exceptionally  well,  adding  both  short  range  interactions  as  well  as  a  long 
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range  Kac  potential,  which  induces  a  liquid-vapor  type  phase  transition  in  such  systems. 
He  was  then  able  to  show  that  this  transition  did  not  show  peculiar  behavior  in  one 
dimension. 

Mr.  Chelst  has  already  published  one  paper  on  his  results  (copy  attached).  Other 
papers  resulting  from  the  thesis  are  now  in  preparation.  A  copy  of  Mr.  Chelst’s  thesis  and 
publication  will  be  sent  separately.  At  the  present  time  Dr.  Chelst  has  a  teaching  position 
at  Rutgers  University. 

Mr.  Kevin  Rosema  was  supported  by  this  grant  for  the  period  of  one  year.  He  is 
currently  continuing  his  graduate  studies  and  expect  to  finish  at  the  end  of  this  year. 
His  work  on  the  grant  involved  fluctuations  in  Coulombic  type  systems.  In  particular  he 
developed  several  computer  program  to  obtain  systems  whose  fluctuations  are  subnormal, 
i.e.  grow  slower  than  the  volume.  This  is  a  characteristic  of  Coulomb  systems  and  possibly 
other  related  systems.  This  work  has  not  been  published. 

Appendix  A:  Three  copies  of  the  manuscript:  Absence  of  Phase  Transitions  in  Modified 
Two-Component  Plasmas:  The  Analytic  Theory  of  Continued  Fractions  in  Statistical  Me¬ 
chanics,  by  Dov  Chelst. 

Appendix  B:  Modified  Two  Component  Plasmas  and  Generalizations  of  Schwarz’s  Lemma, 
by  Dov  Chelst,  Ph.D.  Thesis 
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Contemporary  Mathematics 
Volume  236.  1999 


Absence  of  Phase  Transitions  in  Modified  Two-Component 
Plasmas:  The  Analytic  Theory  of  Continued  Fractions  in 
Statistical  Mechanics 

Dov  Chelst 


Abstract.  In  1961,  A.  Lenard[9]  and  S.  Prager[ll]  independently  solved  the 
one-dimensional  two-component  plasma.  Their  solutions  employed  an  implicit 
equation  involving  a  continued  fraction.  In  addition,  Lenard  used  the  analytic 
theory  of  continued  fractions  to  prove  the  existence  of  an  implicit  solution. 
However,  he  did  not  realize  that  this  theory  could  also  be  used  to  prove  the 
analyticity  of  this  solution  with  the  aid  of  an  implicit  function  theorem. 

We  have  extended  Lenard ’s  analysis  in  [3]  to  include  a  class  of  systems 
which  we  call  modified  two-component  plasmas.  Any  such  system  can 
still  be  described  in  terms  of  an  implicit  continued  fraction  equation.  In  this 
paper,  we  intend  to  show  that  such  a  solution  is  analytic.  Thermodynami¬ 
cally,  this  implies  the  non-existence  of  pressure-dependent  phase  transitions, 
i.e.  transitions  between  two  phases  that  occur  as  pressure  P  and  density  p 
vary  while  the  inverse  temperature  /3  remains  fixed. 


1.  Introduction:  Statistical  Mechanics,  Phase  Transitions  and 
One-Dimensional  Systems 

It  is  not  our  intent  to  develop  new  theorems  regarding  the  convergence  or 
analyticity  of  continued  fractions.  Rather,  we  intend  to  apply  well-known  theorems 
(Theorems  6  and  7)  to  tackle  a  problem  in  statistical  mechanics.  Thus,  this  paper 
will  contain  no  new  results  for  continued  fractions,  just  a  new  application.  However, 
we  believe  this  application  to  be  quite  significant. 

What  is  statistical  mechanics?  Statistical  mechanics  studies  systems  containing 
a  large  number  of  particles.  These  particles  obey  a  set  of  microscopic  dynamical 
laws.  Starting  from  microscopic  particle  interactions,  such  as  pair  potentials,  one 
seeks  to  derive  expressions  for  thermodynamic  macroscopic  quantities  such  as  the 
pressure  P,  as  functions  of  the  inverse  temperature  (3  and  particle  density  p.  Within 
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this  theory,  these  thermodynamic  quantities  correspond  to  either  averages  or  sums 
over  a  tremendous  number  (e.g.  10^°)  of  microscopic  variables.  The  large  number 
of  particles  is  responsible  for  the  statistical  nature  of  the  theory. 

While  developing  a  thermodynamic  portrait  of  a  given  system,  one  naturally 
asks  if  it  exhibits  a  phase  transition.  A  collection  of  particles  can  behave  in  different 
ways:  like  individual  molecules  (a  gaseous  state),  like  a  loosely  bound  conglomerate 
(a  liquid  state),  or  like  a  more  rigidly  bound  molecular  lattice  (a  solid  state).  The 
question  yet  remains:  for  a  given  system,  do  clear  demarcations  exist  between  these 
different  states?  Do  they  represent  distinct  phases  of  matter? 

A  phase  transition  is  characterized  by  an  abrupt  change.  For  example,  when 
we  increase  the  pressure  on  a  gas,  while  maintaining  its  temperature,  its  molecules 
may  condense  to  form  a  liquid.  Thus,  at  a  specific  pressure,  a  sudden  change  in 
density  would  occur.  Up  to  this  point,  the  pressure  and  density  would  continuously 
increase  upon  the  gas. 

To  describe  this  behavior  mathematically,  we  have  a  fundamental  thermody¬ 
namic  relation  which  describes  the  system  as  a  function  of  a  complete  set  of 
thermodynamic  variables.  One  example  of  such  a  relation  is  the  average  entropy 
per  particle  s(e,p)  which  depends  upon  the  mean  energy  per  particle  e  and  the 
particle  density  p  of  a  system.  Another  formulation,  the  chemical  potential  (or 
Gibbs  potential)  p{P,l3),  which  gives  rise  to  this  article’s  continued  fraction,  de¬ 
pends  explicitly  upon  a  system’s  pressure  P  and  inverse  temperature  /?.  All  other 
thermodynamic  quantities  are  calculated  as  partial  derivatives  of  this  fundamental 
relation.  For  example,  the  density  p  of  a  system  for  some  P  and  (3  is  related  to  p. 

by 


(1) 


1  _  dp{P,(3) 
p  dP 


l/pl  l/p2  bp 

A  phase  transition  can  be  defined  as  a  discontinuity  in  one  of  p  s  derivatives. 
For  example,  in  the  above  diagram,  a  liquid-vapor  phase  transition  of  van  der 
Waals  type  can  be  described  by  the  non-differentiability  of  a  pressure-density  curve 
(isotherm).  Below  a  critical  temperature,  there  is  a  certain  pressure  which  cor¬ 
responds  to  two  distinct  densities  pi  and  p2,  representing  matter  in  two  different 
phases.  For  any  intermediate  value  of  the  density,  each  phase  occupies  a  propor¬ 
tional  fraction  of  the  system’s  total  volume.  At  either  end  of  this  interval  along  the 
isotherm,  the  curve  has  a  discontinuity  in  its  first  derivative,  i.e.  a  first-order  phase 
transition. 

While  one  can  often  prove  that  a  physical  system  exhibits  thermodynamic 
behavior,  one  can  rarely  calculate  its  fundamental  relation  exactly  based  solely 
upon  the  knowledge  of  its  microscopic  interactions.  Thus,  to  obtain  exact  results, 
one  needs  to  examine  idealized  models.  One  such  idealization  involves  restricting 
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the  analysis  to  one  rather  than  three  dimensions.  In  one  dimension,  a  number  of 
models  have  exact  solutions. 

Rigorous  proofs  for  the  non-existence  of  phase  transitions  are  another  rare 
commodity  in  statistical  mechanics.  These  theorems,  when  they  exist,  focus  on 
showing  that  some  fundamental  relation  describing  a  given  system  is  analytic.  In 
one-dimension,  the  most  famous  non-existence  proof  certainly  follows  this  line  of 
reasoning.  In  1950,  Van  Hove [5]  showed  that  one-dimensional  systems  with  only 
finite- neighbor  interactions  exhibit  no  phase  transitions. 

Plasmas,  which  consist  of  charged  particles,  interact  via  a  long-range  Coulomb 
interaction  and  are  not  covered  by  Van  Hove’s  theorem.  Nevertheless,  in  1962, 
Lenard  and  S.F.  Edwards[4]  showed,  by  analyzing  the  eigenvalues  of  a  differential 
equation  of  Mathieu  type,  that  the  one-dimensional  two-component  plasma  can 
exhibit  no  phase  transitions.^  This  was  proven  in  a  way  that  avoided  the  continued 
fraction  which  featured  so  prominently  in  Lenard’s  original  thermodynamic  solution 
of  the  system. 

We  can  readily  generalize  the  original  thermodynamic  argument  of  [9]  to  include 
modified  two-component  plasmas.^  One  simply  notes  that  Lenard  and  Prager’s 
original  argument  uses  a  Laplace  transform  technique  reminiscent  of  the  solutions 
of  Tonks  and  Takahashi  for  systems  with  only  hard-cores  or  with  only  nearest- 
neighbor  potentials  respectively.^  Then,  one  replaces  Lenard’s  Laplace  transform 
of  the  identity  with  a  more  complicated  Laplace  transform.  However,  it  is  not  clear 
how  to  generalize  Lenard’s  second  argument,  that  proved  the  non-existence  of  phase 
transitions,  in  so  broad  a  fashion.  Luckily,  this  is  not  necessary. 

2.  Modified  Tavo- Component  Plasmas 

When  expressing  the  thermodynamics  of  this  model,  we  fix  (3  and  the  unit 
charge  a.  Thus,  in  order  to  focus  upon  the  crux  of  the  argument  and  eliminate 
distracting  notation,  we  set  a  =  1  and  p  =  1.  We  therefore  intend  to  suppress  all 
dependence  upon  these  variables,  totally  ignoring  them  in  all  but  two  places:  at 
the  ends  of  Sections  4.3  and  4.4. 

-  -  +-  +  +  ++  -  ++  -  -- 

I - • — • - • — • — • - • - • - • - • - - • - • - •-> 

0 

Modified  two-component  plasmas  contain  an  equal  number  of  positively  and 
negatively  charged  particles  with  unit  charges  -hi  and  --1  respectively.  Each  particle 
resides  on  the  positive  half  line,  and  is  specified  by  a  pair  of  coordinates  (ai^Xi) 
describing  its  charge  and  position.  Pairs  of  particles  interact  via  a  one-dimensional 
Coulomb  potential 

When  considered  in  three-dimensional  space,  we  view  these  particles  as  parallel 
“charged  sheets”  that  lie  perpendicularly  to  the  x-axis  and  each  of  which  has  a 
“charge  density” 

^In  fact,  he  proved  it  for  more  general  multi-component  charged  systems  satisfying  certain 
conditions. 

^The  entire  argument  can  be  found  in  [3]. 

^See  [13,  pp. 48-53].  A  translation  of  Takahashi’s  original  article  can  be  found  in  [10,  p.25]. 
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In  addition,  nearest-neighbor  particles,  separated  by  a  distance  x,  interact  via 
a  potential  ipix).  This  potential  may  depend  upon  the  relative  charges  of  the  two 
particles.  In  this  case,  we  refer  to  two  potentials,  ip  same  and  ipopp,  as  the  interactions 
between  neighboring  pairs  of  ions  with  the  same  charges  and  with  opposing  charges 

respectively.  ... 

The  total  potential  energy  H  of  this  system  in  a  given  configuration  is  just 

the  sum  of  all  the  pair  potentials.  Using  H,  for  a  fixed  pressure  P ,  we  employ  a 
statistical  ensemble  especially  suited  to  one-dimensional  calculations,  the  isobaric- 
isothermal  ensemble^,  to  determine  the  fundamental  relation  ii{P)  exactly.  That  is 
all  we  will  say  about  the  basic  description  of  modified  two-component  plasmas  and 
their  statistical  mechanical  calculations.  A  full  derivation  of  n{P)  can  be  found  in 
[3].  We  will  now  proceed  to  the  fruits  of  this  analysis. 

3.  Thermodynamic  Results:  The  Real  and  the  Complex 

The  chemical  potential  p  of  a  modified  one-dimensional  two-component  plasma 
is  related  to  the  implicit  solution  of  the  equation 

(2)  QiPz)  =  l, 

for  positive  P.  Specifically,  the  implicit  function  z*{P)  defined  to  be  the  solution 
of  smallest  modulus  of  (2)  for  fixed  P  >  0,  is  related  to  /x  by 

(3)  P(P)  =  e2^(^). 

Through  this  relation,  we  will  show  that  ji(P)  is  an  analytic  functioii  of  P. 

Before  we  can  examine  we  need  more  information  about  Q.  Q  can  be 
described  in  two  ways: 

1.  Power  Series:  Initially,  for  fixed  P  >  0,  we  define  Q  as  a  power  series  in  z 
about  0, 

oo 

(4)  = 

n=:l 

Each  coefficient  Qn{P)  is  the  result  of  a  statistical  calculation  involving 
a  modified  two-component  plasma  with  2N  particles  using  the  isobaric- 
isothermal  ensemble.  As  a  consequence,  each  coefficient  is  a  Laplace  trans¬ 
form  of  a  nonnegative  function^,  whose  integral  converges  for  all  positive 
pressures  P  and  hence  for  all  complex  P  with  positive  real  part.  Our  main 
concern  is  that  it  is  positive  for  all  positive  P. 

2.  Continued  Fraction:  Later  in  [3],  we  show  that  after  imposing  a  few 
conditions  and  employing  a  recursive  analysis,  a  modified  two-component 
plasma’s  Q  can  be  described  by  a  continued  fraction  whose  partial  numera¬ 
tors  and  denominators  are  functions  of  P  and  z. 

The  first  description  shows  us  that  when  P  >  0,  and  z  >  0  is  within  the  radius 
of  convergence  about  0,  Q(P,  z)  and  ^  are  positive.  The  second  description  allows 
us  to  fully  utilize  complex  analytic  machinery.  With  it,  we  can  show  that  Q  is  a 
meromorphic  function  of  both  P  and  z.  We  will  see  that  Q’s  continued  fraction  is 

^See  Percus[8]  for  a  discussion  of  this  ensemble. 

5q^(P)  corresponds  to  a  portion  of  the  isobaric-isothermal  partition  function  and  is  thus  a 
Laplace  transform  of  some  portion  Zm{L)  >  0  of  the  canonical  partition  function. 


CONTINUED  FRACTIONS  AND  MODIFIED  TCP’S 


141 


well-defined  when  P  is  restricted  to  some  appropriate  open  set  containing  the 
positive  reals  R+,  while  z  is  allowed  to  range  over  all  C. 

The  continued  fraction  describing  Q  can  be  written  in  the  standard  form 


(5)  (3(P,z)  =  K(^j. 

When  Ip  is  charge-dependent,  the  partial  numerators  and  denominators  are  quite 
complicated: 


(6) 


®(-^  +  1)> 

Cln  ~  —Z  Q{P  +  71^), 


bn  —  VoppiP  4"  n^)  —  z  c,(P  -f  (n  -t- 1)^). 
Each  Tj  is  the  reciprocal  of  a  Laplace  transform: 

1 


=  £(e-^> 


=  C{e  and 


Vsameis)  ''  r]opp{s)  rj,,g{s) 

E  and  ©  depend  upon  rjopp  and  risame  through  the  relations 


-I- ' 


Vsamei^)  ^opp{^) 


(7) 


P)/  s  _  /  Voppis) 

^  ’  Wameis) 


,  and  E  = 


Voppis) 


(1-6(5)). 


When  tp  does  not  depend  upon  the  neighboring  particles’  charges,  this  descrip¬ 
tion  simplifies  drastically.  V'opp,  '>Psame  and  ipavg  all  become  equal,  while  0  =  1  and 
E  =  0.  Hence, 


(8) 


'“1  - 
On  = 

bo  =  0,  and 
,b„  =  r?(P-f  n^). 


The  fact  that  the  description  of  these  continued  fractions  hinges  upon  the 
reciprocal  of  a  Laplace  transform,  coupled  with  the  desire  to  have  the  continued 
fraction  description  remain  valid  for  all  positive  P,  motivates  us  to  stipulate  the 
following  condition  for  ip. 


Condition  1.  The  Laplace  transform  C{e~'^){P)  determined  by  any  nearest 
neighbor  potential  ip  describing  a  modified  two-component  plasma  will  converge  for 
all  positive  P  and  hence  in  the  complex  right  half  plane  Hq  =  {a;  :  8?(z)  >  0}. 

This  condition  is  not  specific  to  modified  plasmas;  it  is  implicit  in  the  standard 
thermodynamic  description  of  pure  nearest  neighbor  systems  solved  by  Takahashi[10, 
p.25].  Thus,  Condition  1  is  imposed,  implicitly  or  explicitly,  on  many  systems  con¬ 
taining  nearest-neighbor  potentials. 

Of  course,  we  must  now  convince  a  reader  that  the  reciprocals  of  Laplace  trans¬ 
forms,  the  77’s,  are  well-defined  for  appropriate  values  of  P.  In  other  words,  we  need 
to  show  that  the  Laplace  transforms  are  non-zero.  Only  then  can  we  show  that  the 
continued  fractions  given  by  (6)  and  (8)  converge.  We  will  explore  the  properties 
of  Laplace  transforms  which  make  each  r?  well-defined  in  Section  4.3.2. 
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4.  The  Analyticity  of  n{P) 

4.1.  The  Main  Result.  Now,  we  will  state  our  main  theorem  and  prove  it 
assuming  Lemma  3  that  we  will  prove  in  Section  4.4. 

Theorem  2.  Let  a  nearest-neighbor  potential  ip  satisfy  Condition  1.  When  ip 
is  charge- dependent,  let  it  satisfy  Condition  8.  Let  Q{P,z)  he  a  function  that  can  be 
described  both  as  a  continued  fraction  (6)  and  as  a  power  series  (4)  for  fixed  P>Q 
in  z  about  0  with  positive  coefficients.  Then  the  implicit  function  ffiP)  defined  via 
(2)  and  (3)  is  analytic  for  P  >  0. 

4.2,  Proof  of  Theorem  2.  To  prove  Theorem  2,  we  must  traverse  a  number 
of  steps. 

Proof.  1.  First,  we  state  the  lemma; 

Lemma  3.  IfQ{P,z)  satisfies  the  hypotheses  of  Theorem  2,  then  Q  is 
a  separately  meromorphic  function  of  P  and  z  for  all  {P,  z)  &  x  C- 
for  some  appropriate  open  subset  of  C  containing  M+ . 

IVhile  we  leave  the  proof  of  Lemma  3  to  Section  4.4,  we  would  like  to 
describe  its  conclusion  and  give  an  idea  of  its  proof  here.  Specifically,  we 
show  that  if  we  fix  one  variable,  and  examine  a  neighborhood  about  any 
value  of  the  other  variable,  there  is  an  appropriate  fcth  remainder  of  the 
continued  fraction,  K(on+A:/f>n+fc)!  which  is  analytic  in  that  neighborhood. 
Thus,  the  full  continued  fraction  is  a  rational  function  of  the  first  k  partial 
numerators  and  denominators  and  the  remainder  term.  Since  these  are  all 
analytic  functions  of  the  free  variable,  Q  must  be  a  meromorphic  function 
of  the  free  variable. 

This  argument  utilizes  two  theorems  from  the  analytic  theory  of  contin¬ 
ued  fractions  quoted  in  Section  4.4  as  Theorems  6  and  7.  At  the  same  time, 
it  relies  heavily  upon  the  asymptotic  behavior  of  Laplace  transforms  de¬ 
scribed  in  Section  4.3.  It  is  not  suprising  that  for  a  charge-dependent  system, 
we  need  an  extra  condition,  given  in  Section  4.4  as  Condition  8,  to  apply 
these  continued  fraction  theorems.  A  charge-independent  system  does  not 
require  this  additional  condition;  yet,  when  ip  assumes  the  proper  form,  it 
automatically  satisfies  Condition  8. 

2.  An  extension  of  Hartogs’  Theorem  to  meromorphic  functions  due  to  W. 
Rothstein[12]  states  that  a  separately  meromorphic  function,  e.g.  Q,  in  an 
open  complex  domain,  e.g.  x  C,  is  actually  jointly  meromorphic.  Thus, 
Q  is  jointly  analytic  in  a  neigborhood  of  any  solution  to  (2). 

3.  Moreover,  we  note  that  the  power  series  description  (4)  of  Q  implies  that 
^  is  non-zero  about  any  solution  to  (2)  provided  that  both  coordinates 
are  positive  and  z  lies  in  the  radius  of  convergence  of  Q.  Together  with  the 
meromorphicity  of  Q,  this  assures  us  that  for  P  >  0,  Q's  first  pole  away  from 
zero  occurs  on  the  positive  real  axis.  It  also  guarantees  that,  for  positive  P, 
a  unique  positive  solution  z{P)  exists  within  this  radius  of  convergence 
and  that  this  solution  has  the  smallest  modulus  of  all  solutions  to  (2).  Thus, 
this  solution  is  actually  z*(P). 

4.  We  may  now  employ  the  complex  analytic  implicit  function  theorem[6, 
1.B.6]  to  prove  that  z*(P)  is  an  analytic  function  of  P.  Let  us  briefly  review 
the  implicit  function  theorem. 
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Theorem  4  (Implicit  Function  Theorem).  Let  f  :  U  C  be  a  jointly 
analytic  function  in  some  open  domain  U  C  In  addition,  let  f{x,  z)  ~  0 
for  some  (f,  z)  e  U ,  and  |^(f,  z)  0.  Then  there  exists  a  neighborhood  U' 
about  X  and  an  analytic  function  g  :  — >  C  such  that  f{x,g{x))  =  0  and 

(x,g{x))  e  U  for  all  x  e  U' . 

This  theorem  applies  to  our  situation  when  we  replace  /  by  Q  “  1  •  Since 
/  is  meromorphic,  it  is  automatically  analytic  in  a  neighborhood  of  any  point 
that  is  not  a  pole  of  /  which  includes  any  point  where  /  =  0.  Thus,  by  the 
implicit  function  theorem,  z*(P)  is  analytic  in  an  appropriate  neighborhood 
of  P. 

5.  Since  P  >  0  was  chosen  arbitrarily,  z*(P)  must  be  analytic  for  all  positive 
P.  The  fact  that  the  solution  of  smallest  modulus  z*{P)  is  always  real 
and  positive  makes  this  global  result  possible.  Normally,  given  two  analytic 
implicit  solutions  Zi{P)  and  Z2{P)  to  a  single  equation,  a  situation  might 
arise  in  which  z\  would  have  the  smaller  modulus  for  some  values  of  P  while 
Z2  would  have  the  smaller  modulus  for  other  values  of  P.  The  fact  that  the 
solution  ^*(P)  of  smallest  modulus  is  always  positive  implies  that  for  such 
a  switch  to  take  place  Zi{P)  and  Z2[P)  must  coincide  for  some  value  of  P. 

By  the  assumptions  of  the  implicit  function  theorem,  namely  that  ^  7^  0, 
this  cannot  occur. 

6.  The  fact  that  z*(P)  is  analytic  and  positive  allows  us  to  conclude  that 
li[P)  =  I  ln(z*(P))  is  also  analytic  in  P. 

□ 

4.3.  Three  Features  of  Laplace  Transforms.  As  Q  is  defined  in  terms  of 
Laplace  transforms,  we  should  review  their  relevant  features.  More  specifically,  we 
list  three  features  of  Laplace  transforms.  The  first  is  satisfied  for  any  function  with 
an  absolutely  convergent  Laplace  integral.  The  second  is  specific  to  transforms 
of  nonnegative  functions  (that  are  not  identically  zero),  which  is  certainly  true  in 
our  case  as  e~'^  is  never  negative.  The  third  holds  true  only  for  an  even  more 
specific  situation;  this  is  meant  to  serve  as  an  example  of  a  Laplace  transform’s 
leading-order  asymptotics. 

4.3.1.  Limiting  Behavior.  If  the  Laplace  transform  of  a  function  /,  C{f){s)  = 

/o^  f{t)e~^^dt  converges  absolutely  for  some  s,  it  is  analytic  and  defined  for  all  P  in 
some  open  right  half  plane  =  {s  :  9f^(s)  >  w}.  In  addition,  >^(/)(*5)  = 

O.  In  our  situation,  by  the  assumptions  of  Condition  1,  C{e~'^){s)  is  defined  in  Ho 
and  for  s  €  C  Ho,  limT^^oo  \v{s  +  n^)|  =  00. 

4.3.2.  Transforms  of  Nonnegative  Functions  and  Their  Zeros.  For  a  nonneg¬ 
ative  function  /  that  converges  in  Hq,  this  limiting  behavior  is  monotonic  for  real 

P.  In  addition,  for  complex  s,  \r]{s)\  >  r){^{s)).  Moreover,  if  /  is  nonnegative  and 
not  identically  0,  we  can  be  sure  that  rj  is  well-defined  in  some  neighborhood  O 
containing  IR+,  After  all,  T(/)(s)  is  anal3d:ic  in  Ho  and  strictly  positive  on  R_i_. 

Unfortunately,  this  does  not  suffice.  It  is  necessary  to  find  a  complex  neighbor¬ 
hood  Oe,  containing  M+,  in  which  C{f){s)  and  C{f){s  +  n^)  are  nonzero  for  all  n. 
We  will  prove  that  there  exists  a  neighborhood  containing  E+,  on  which  £(/) 
is  nonzero,  and  which  has  the  property  that  s  e  implies  s-\-5  £  all  J  >  0. 

In  other  words: 
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Lemma  5.  Let  C{f)  be  a  Laplace  transform  of  a  nonnegative  function  f  that 
converges  for  all  P  G  Hq.  If  f  is  not  identically  zero,  then  the  zeros  of  /l[f)  are 
bounded  away  from  M+  in  any  half-plane  for  x  >  0. 


Proof.  We  will  show  more  than  the  lemma  requires.  We  will  show  that  the 
zeros  of  the  real  part  of  C{f)  are  also  uniformly  bounded  away  from  the  positive 
real  axis  in  Ux  for  any  x  >  0.  Given  the  original  function  /  and  a  positive  real 
value  X,  pick  some  e  <  |.  Now,  by  simple  convergence  arguments,  one  can  choose 

a  constant  K  so  that  f{t)e~^^dt  >  0  and 


poo 

/  f{t)e-^^dt  <  e  /  fit)e-^*dt. 

Jk  Jo 


Then,  for  all  si  >  x, 

pOO  /’OO 

/  fit)e-^^*dt  <  /  fit)e-^^dt  < 

Jk  Jk 


<  e 


r  f{t)e--^dt 
Jo 

[  f{t)e-^^^dt 
Jo 


Now,  simply  choose  it;  >  0  so  that  wK  <  f .  For  all  real  values  |52|  <  w  and 
0  <t  <  K,  cos  S2t  >  ^  and  for  complex  s  =  si  +  is2  in  the  complex  strip  Sx  =  {s  : 

R{s)  >  X,  |^(s)|  <  w}, 

1  pK  poo 

sR{Z:(/)(5)}=  /  fit)coss2te-^^^dt>-  fit)e-^^^dt-  dt 

Jo  ^  Jo  ^K 

.K 

I  /  /(<)e 

Jo 


>  0. 


The  proof  is  complete  when  we  note  that  since  Sj,  contains  no  zeros  of  T(/),  the 
Laplace  transform’s  zeros  must  all  remain  at  least  a  distance  «;  >  0  away  from  R+ 
in  H.. 


Before  we  continue,  we  simply  construct  an  open  set  using  the  conclusion  of 
the  previous  lemma.  Given  a  suitable  nonnegative  function  /  and  a  positive  e  <  2 ! 
define  the  “zero-less”  open  set  =  Ux>o  Thus,  because  C{f)  has  no  zeros  and 
is  analytic  in  we  are  assured  that  rj  is  well-defined  and  analytic  there  as  well. 
Finally,  we  note  that  this  lemma  completes  a  pair  of  inequalites.  For  any 

S  € 

(9)  |£(/)(s)|  <  T(/)TO)  <  \f^\C{f){s)\. 

2  ^ 

4.3.3.  Leading  Order  Asymptotics  of  C{f)  -  Watson's  Lemma.  While  it  is  not 
completely  necessary,  the  inequality  (9)  allows  us  to  consider  the  complex  asymp¬ 
totic  behavior  of  C{f){s)  as  5R(s)  oo.  While  the  limit  of  £(/)  is  certainly  an 
asymptotic  feature,  in  a  charge- dependent  modified  plasma  system,  where  two 
Laplace  transforms  play  a  role,  we  need  a  clearer  understanding  of  Liffs  asymp¬ 
totic  behavior. 

Toward  this  end,  we  make  a  large  assumption  about  / .  Every  nonnegative  func¬ 
tion  has  a  set  /"^((O,  oo))  €  M-f  of  points  where  /  is  positive.  For  a  function  /  that 
has  a  right-handed  limit  near  the  smallest  point  p  in  this  set’s  closure,  £(/)(s)’s 
asymptotic  behavior  for  large  5?(s)  depends  very  strongly  upon  this  limit.  Watson  s 
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Lemma[2,  6.4]  covers  a  number  of  possible  limiting  behaviors.  While  this  does  not 
exhaust  all  possible  functions,  nor  does  it  include  all  physically  meaningful  func¬ 
tions^,  we  will  restrict  ourselves  to  those  potentials  covered  by  Watson’s  Lemma. 
An  industrious  reader  will  see  that  those  generalizations  of  Watson’s  Lemma  that 
are  also  found  in  [2]  can  be  treated  similarly. 

Watson’s  Lemma  states  that  if  a  function  /  behaves  as  At^  to  leading  order  as 
t  ^  0“^,  for  some  constants  A  and  a  >  0  and  if  £(/)  converges  in  some  then 

the  Laplace  transform  C{f){s)  decays  as  ^{s)  oo,  to  leading  order,  like 
(r  is  the  Gamma  function  which  for  positive  integers  n  satisfies  r(n  +  1)  =  n\  and 
using  Euler’s  integral  can  be  written  as  r(x)  =  dt. )This  is  certainly  true 

for  real  s  and  the  statement  for  complex  s  follows  from  (9). 

In  statistical  mechanics,  we  frequently  encounter  potentials  '0  that  contain  hard 
core  exclusions.  In  terms  of  /  =  a  hard  core  specifies  that  f{t)  —  0  for  all 

t  <  b  for  some  hard  core  diameter  b.  If,  in  addition,  as  ^  ^  /  behaves  as 

A{t  -  b)^  for  some  constants  A  and  a  >  0,  its  Laplace  transform  C{f){s)  decays, 

to  leading  order,  as  - ■  This  follows  from  Watson’s  Lemma  and  an 

elementary  feature  of  Laplace  transforms  that  can  be  found  in  any  undergraduate 
text  on  the  subject.  The  reciprocal  of  C{f)  in  either  of  these  two  situations,  6  =  0 
or  6  /  0, diverges  to  infinity  in  an  appropriate  manner. 

We  note  that  while  Watson’s  Lemma  also  allows  for  -1  <  a  <  0,  a  physical  ar¬ 
gument  disallows  these  exponential  values.  After  all,  since  /  =  our  exponent 
a  depends  linearly  upon  the  inverse  temperature  (3.  Thus,  while  for  some  values 
of  /?,  0  >  a  =  a'/?  would  be  greater  than  -1,  there  would  certainly  be  values  of  /? 
for  which  this  would  be  false  and  thus  the  Laplace  transform  would  not  even  exist. 
This  contradicts  Condition  1  that  we  imposed  upon  our  nearest-neighbor  potential 
0  in  Section  3. 

This  section  on  Laplace  transforms  provides  all  that  we  need  in  this  paper. 
It  furnishes  analyticity  criteria.  It  gives  a  complex  open  set  in  which  all  77’s  are 
well-defined,  i.e.  just  the  finite  intersection  of  open  sets  for  each  separate  r/.  It  also 
provides  the  asymptotic  information  regarding  our  ry’s  that  allows  us  to  apply  our 
continued  fractions  theorems. 

4.4.  The  Meromorphicity  of  Q. 

Proof  of  Lemma  3.  As  stated  in  the  proof  of  Theorem  2,  our  aim  is  to  show 
that  some  remainder  of  the  continued  fraction  Q  is  separately  analytic  in  x  C; 
i.e.  analytic  in  a  single  variable  while  the  other  remains  fixed.  For  this  purpose,  we 
employ  two  theorems  found  in  chapter  4  of  the  book  by  Thron  and  Jones  [7,  Thms. 
4.35  and  4.54].  Note  that  refers  to  the  nth  convergent  of  the  continued  fraction 
whose  value  we  can  obtain  by  substituting  0  for  Un+i- 

Theorem  6  (Pringsheim) .  //  |bn|  >  |an|  +  1  for  all  n,  then 

i.  K  j  converges,  and 

ih  Ifni  <  1  for  all  n. 


^Trye 


from  a  potential 
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Theorem  7.  If: 

i.  On (2:),  bn{x)  are  analytic  in  some  domain  D,  and 

ii.  3^1,6  €  C,  s.t.  for  all  x  e  D  and  n  €  Z+,  f„(a;)  ^  {^i,6)0o}- 

Then  K  *5  analytic  in  D. 

We  first  note  that  all  the  partial  numerators  and  denominators  described  in  (6) 
and  (8)  are  analytic  functions  of  P  and  z.  In  O^,  both  rjopp  and  r/same  nonzero 
and  hence  their  ratio  is  an  analytic  function.  Thus,  0  and  S  are  both  analytic 
functions  of  P.  Analyticity  in  z  is  even  more  obvious. 

Now,  consider  e(P)  and  E{P)  defined  in  (7).  By  the  asymptotic  behavior  of 
the  two  7?’s,  according  to  Watson’s  Lemma, 


|0(s)| 


2(Qopp  +  1)  2bopp50 
C.^0  e  - 


2(asame  +  l)  _26 


=  Csl^  Otopp- 


e  )  g2  (feopp  ^same  and 


^0 


DSr 


ame  50 

^(^opp 

same  )So 


for  some  constants  C  and  D  when  5  is  in  and  5o  —  3^(5). 

Since  E{s)  can  never  grow  asymptotically  larger  than  0(s),  if  we  wish  to  prove 
that  671  is  eventually  bigger  than  an  we  must  focus  on  its  other  term,  T]opp{P  +  n.  ) 
which  grows  like  with  sq  =  ^{P)  +  n^.  For  this  to  grow  more 

quickly  than  its  competing  partial  numerator,  we  impose  the  following  condition. 


Condition  8.  Let  ^  be  a  nearest-neighbor  interaction  within  a  modified  two- 
component  plasma.  Moreover,  let  /(t)  =  e  have  the  form  f{t)  =0  for  t  <  b 
and  f{t)  as  f  6+  for  some  nonnegative  values  b  and  a.  Then,  we 

require  that 

1.  bopp  ^  2650,7167  and 
ii.  if  bppp  2650,716,  that  oc^pp  ^  2oi5077ie' 


If  this  condition  holds,  we  are  assured  that  eventually  bp  will  grow  asymptoti¬ 
cally  larger  than  an  and  the  conditions  of  Pringsheim’s  lemma  (Lemma  6)  will  be 
satisfied.  After  all,  whether  we  fix  P  or  z,  the  real  part  of  P  -f  n^,  that  replaces  s 
from  the  previous  discussion,  goes  off  to  infinity  and  our  asymptotic  analysis  holds. 
The  second  lemma  will  automatically  follow  if  we  choose  any  two  complex  numbers 
and  $2  with  modulus  greater  than  1.  Consequently,  Q{P,z)  is  meromorphic  in 


Remark  9.  We  would  like  to  point  out  that  charge- independent  systems  do 
not  require  the  additional  condition  stated  above.  If  one  examines  their  partial 
numerators  in  (8),  one  notices  that  they  remain  fixed  as  n  increases.  Conversely, 
the  partial  numerators  grow  asymptotically  large  regardless  of  the  special  form  of 
e~^.  Thus,  both  lemmas  provide  their  conclusions  effortlessly. 

Remark  10.  In  addition,  while  we  have  not  ruled  out  temperature-dependent 
phase  transitions,  it  is  possible  that  a  similar  analysis  will  preclude  their  existence 
as  well.  Paying  attention  to  the  inverse  temperature  (5  merely  involves  replacing 
pjpn'^hy  /?(P-f  n^)  in  77’s  arguments  and  switching  the  transformed  function  from 
to  The  dependence  upon  (3  that  arises  from  multiplying  the  arguments 

of  r?  by  /3  can  be  treated  in  exactly  the  same  manner  that  we  treated  P  above.  Only 
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the  (3  which  figures  in  the  definition  of  the  transformed  function  itself  e  could 
prove  problematic. 

In  certain  simple  instances,  e.g.,  when  ^  =  0  or  when  ^  is  a  pure  hard  core 
exclusion,  temperature-dependent  phase  transitions  are  also  absent.  In  those  situ¬ 
ations,  C{e~'^)  depends  on  (3  only  through  its  argument  (3{P  -h  n^).  Thus,  in  these 
instances,  the  line  of  reasoning  described  in  this  paper  applies  to  jd  as  well. 

5.  Conclusion 

We  have  accomplished  our  goal  of  utilizing  the  theory  of  continued  fraction  to 
prove  a  fact  regarding  modified  two-component  plasmas.  Of  course,  some  questions 
remain  unresolved.  Did  we  use  this  theory  effectively?  The  theorems  employed 
provided  sufficient  conditions  preventing  the  existence  of  phase  transitions  in  these 
systems.  But,  how  necessary  is  the  additional  condition  imposed  upon  plasmas  with 
charge-dependent  nearest  neighbor  interactions.  What  happens  if  bopp  >  2bsame^ 

Continued  fractions  may  yet  hold  the  answer  to  this  question.  Continued  frac¬ 
tions  continue  to  persist  in  many  models  related  to  the  two-component  plasma. 
They  appear  explicitly  in  the  discussion  of  one-dimensional  ion-dipole  systems  in 
[15].  Moreover,  they  are  implicit  in  every  article  that  calls  upon  the  reader  to  calcu¬ 
late  an  eigenvalue  of  a  Mathieu-type  differential  (differential-delay)  equation  with 
periodic  boundary  values.’^  For  these  reasons,  we  believe  that  continued  fractions 
will  continue  to  provide  valuable  information  regarding  this  family  of  statistical 
models. 
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